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[57] ABSTRACT 

An ultrasonic wave is applied to an anisotropic sample 
material in an initial direction and an angle of flux devia- 
tion of the ultrasonic wave front is measured from this 
initial direction. This flux deviation angle is induced by 
the unknown applied load. The flux shift is determined 
between this flux deviation angle and a previously de- 
termined angle of flux deviation of an ultrasonic wave 
applied to a similar anisotropic reference material under 
an initial known load condition. This determined flux 
shift is then compared to a plurality of flux shifts of a 
similarly tested, similar anisotropic reference material 
under a plurality of respective, known load conditions, 
whereby the load applied to the particular anisotropic 
sample material is determined. A related method is 
disclosed for determining the fiber orientation from 
known loads and a determined flux shift. 

14 Claims, 3 Drawing Sheets 
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METHODS OF DETERMINING LOADS AND 
FIBER ORIENTATIONS IN ANISOTROPIC 
NON-CRYSTALLINE MATERIALS USING 

ENERGY FLUX DEVIATION 5 

ORIGIN OF THE INVENTION 

The invention described herein was made by employ- 
ees of the United States Government and may be manu- 
factured and used by or for the Government for govern- 10 
mental purposes without the payment of any royalties 
thereon or therefor. 

BACKGROUND OF THE INVENTION 

1. Technical Field of the Invention 15 

The present invention relates generally to load moni- 
toring and more particularly to a method of determin- 
ing loads and fiber orientations in anisotropic materials 
using energy flux deviation. 

2. Discussion of the Related Art 20 

In isotropic materials, the direction of the energy flux 

(energy per unit time per unit area) of an ultrasonic 
plane wave is always along the same direction as the 
normal to the wave front. In anisotropic materials, how- 
ever, this is true only along symmetry directions. Along 25 
other directions, the energy flux of the wave deviates 
from the intended direction of propagation. This phe- 
nomenon is known as energy flux deviation and is illus- 
trated in FIG. 1. An ultrasonic transducer U is coupled 
to an anisotropic crystalline material AM and directs an 30 
ultrasonic wave through the material as indicated by 
directional arrow N which is normal to the wave front. 
The anisotropic nature of the material causes the energy 
flux to deviate from the normal arrow N by an energy 
flux deviation angle A, resulting in an energy flux vec- 35 
tor V associated with the deviated wave flow W. The 
direction of the energy flux is dependent on the elastic 
coefficients of the material. This effect has been demon- 
strated in many anisotropic crystalline materials. In 
transparent quartz crystals, Schlieren photographs have 40 
been obtained which allow visualization of the ultra- 
sonic waves and the energy flux deviation. 

The energy flux deviation in graphite/epoxy (gr/ep) 
composite materials can be quite large because of their 
high anisotropy. The flux deviation angle has been cal- 45 
culated for unidirectional gr/ep composites as a func- 
tion of both fiber orientation and fiber volume content. 
Experimental measurements have also been made in 
unidirectional composites. It has been further demon- 
strated that changes in composite materials which alter 50 
the elastic properties such as moisture absorption by the 
matrix or fiber degradation can be detected nondestruc- 
tive^ by measurements of the energy flux shift. 

Graphite fiber-reinforced composites such as gra- 
phite/epoxy, graphite/magnesium and gra- 55 
phite/aluminum exhibit very high stiffness-to-weight 
and strength-to-weight ratios, making them excellent 
materials for lightweight aerospace structures. Since 
these structures are intended primarily to carry load, 
the in-situ and non-destructive determination of load 60 
and load induced quantities such as stress is very desir- 
able. The energy flux of stress waves, i.e., waves af- 
fected by a load, propagating through anisotropic crys- 
tals has been shown to deviate from the direction of the 
normal to the plane wave. However, there has been no 65 
known indication that this deviation can be correlated 
with the amount of applied load. This lack of knowl- 
edge is not surprising since crystals are normally not 


used in load bearing applications. In addition, there is no 
known work concerning any energy flux deviation in 
highly anisotropic, non-crystalline materials such as 
graphite fiber-reinforced composites and consequently 
no known work concerning any effect of applied load 
on such deviations. 

OBJECTS 

It is accordingly an object of the present invention to 
determine load and load induced quantities such as 
stress and strain in anisotropic materials using energy 
flux deviations. 

It is another object of the present invention to accom- 
plish the foregoing object non-destructively. 

It is a further object of the present invention to ac- 
complish the foregoing objects in a simple manner. 

Additional objects and advantages of the present 
invention are apparent from the specification and draw- 
ings which follow. 

SUMMARY OF THE INVENTION 

The foregoing and additional objects are obtained by 
a method of determining a load applied to a particular 
anisotropic sample material according to the present 
invention. An ultrasonic wave is applied to the aniso- 
tropic sample material in an initial direction and an 
angle of flux deviation of the ultrasonic wave front is 
measured from this initial direction. This flux deviation 
angle is induced by the unknown applied load. The flux 
shift is determined between this flux deviation angle and 
a previously determined angle of flux deviation of an 
ultrasonic wave applied to a similar anisotropic refer- 
ence material under an initial known load condition. 
This determined flux shift is then compared to a plural- 
ity of flux shifts of a similarly tested, similar anisotropic 
reference material under a plurality of respective, 
known load conditions, whereby the load applied to the 
particular anisotropic sample material is determined. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a prior art device used to show the energy 
flux deviation which is dependent on the elastic coeffici- 
ents of the material; 

FIG. 2 is an apparatus according to the present inven- 
tion for determining loads in anisotropic non-crystalline 
materials using energy flux deviation; 

FIG. 3 is a representative graph correlating the signal 
ratio of two adjacent receiving transducers and the 
applied load; 

FIG. 4 is a graph of the energy flux deviation as a 
function of the fiber orientation angle in a no applied 
stress condition; 

FIG. 5 is a graph of the shift in energy flux deviation 
as a function of the fiber orientation angle due to a 1 
GPa load applied along the X3 axis, i.e., in the fiber 
direction; and 

FIG. 6 is a graph of the shift in energy flux deviation 
as a function of tie fiber orientation angle due to a 0. 1 
GPa load applied along the xi axis, i.e., in the laminate 
stacking direction. 

DETAILED DESCRIPTION 

Referring to FIG. 2, a transmitting ultrasonic trans- 
ducer 10 is acoustically coupled to a non-crystalline, 
highly anisotropic material AM, e.g., a graphite fiber- 
reinforced composite such as graphite/epoxy, gra- 
phite/magnesium or graphite/aluminum. The ultra- 
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sonic transducer 10 propagates an ultrasonic wave 
through the anisotropic material AM which initially has 
a wave front which is normal to directional arrow N. 

As discussed in greater detail below, an applied load LI 
and L2 causes energy flux deviation angles of A°QL and 5 
A°QT. These propagation directions change with the 
respective applied loads. The directional arrows QL 
and QT are normal to the deflected wave fronts. 

The following variables are used throughout the de- 
scription of the invention: 10 

6 angle, measured from the normal vector N, of fiber 
orientation of anisotropic material AM; 

QT vector designating “quasi-transverse” particle dis- 
placement of the propagated ultrasonic wave; 

QL vector designating “quasi-longitudinal” particle 15 
displacement of the propagated ultrasonic wave; 

A°QT angle, measured from the normal vector N, of 
energy flux deviation in QT direction; and 
A°QL angle, measured from the normal vector N, of 
energy flux deviation in QL direction. 20 

An array of receiving transducers respectively la- 
beled Vi, V2, V3, . . . Vivare acoustically coupled to the 
non-crystalline, highly anisotropic material AM oppo- 
site transmitting transducer 10. The receiving transduc- 
ers are arranged such that adjacent transducers can 25 
detect the desired increment of propagation shift, i.e., of 
the deviation of energy flux, of the wave. Such an array 
is disclosed in U.S. Pat. No. 4,499,770, the specification 
of which is hereby incorporated by reference. The 
transmitting and receiving transducers can be any con- 30 
ventional type such as piezoelectric or capacitive sen- 
sors. The receiver can alternatively be a scanned optical 
ultrasonic receiver or any array of such receivers. 

As known, the received ultrasonic wave is converted 
by each receiving transducer into a respective electrical 35 
signal having a voltage proportional to the wave inten- 
sity. In one embodiment, the shift in propagation caused 
by an applied load is monitored as a simple ratio be- 
tween two adjacent receiving transducers, e.g., V2:V3. 

A sample correlation between the applied load and this 40 
signal ratio is shown in FIG. 3. Material specific data 
curves are used to determine the applied load on that 
material with the measured signal ratio. 

In a more specific application, the effects of nonlinear 
elasticity on energy flux deviation in unidirectional 45 
composites such as graphite/epoxy (gr/ep) was consid- 
ered. The angle of the energy flux deviation, AQT or 
AQL, for such a material was shown to be a function of 
applied stress as a result of this elastic nonlinearity. This 
deviation was modeled as detailed below using acousto- 50 
elastic theory and second and third order elastic stiff- 
ness coefficients for T300/5208 gr/ep previously deter- 
mined in NASA Contractor Report 4100, “Ultrasonic 
Characterization of the Nonlinear Elastic Properties of 
Undirectional Graphite/Epoxy Composites”, William 55 
H. Prosser, Oct., 1987. 

Assuming linear elasticity, the components of the 
energy flux vector (E/) are a function of the linear elas- 
tic stiffness coefficients Cijkl and the spatial and time 
derivatives of the displacement vector (u/). The equa- 60 
tion is given by 

where the Einstein summation convention on repeated 
indices is assumed throughout this discussion. The angle 
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of energy flux deviation can then be calculated as the 
angle between the flux vector and the normal to the 
plane wave front. 

To include non-linear elastic effects on the energy 
flux deviation, acoustoelastic theory is used. This the- 
ory predicts an “effective” linear elastic stiffness tensor 
(c *ijki) that is a function of the second and third order 
elastic coefficients and the applied stress (tr,y), as de- 
scribed in a communication to applicants from Dr. 
David M. Barnett of Stanford University. The expres- 
sion for the “effective” stiffness tensor is given by 

c*nlij=k n Uj+cr„j8li ( 2 ) 

where Bn is the Kronecker delta and k n iij is given by 

knlij = C n Uj + Crlijtnr + Cnstfi-is + Cnlpf-ip + Cnliqtjq + 0) 

Cnlijuvtuv + Crlijut^ui^nr + C/isiJut^uv^is + 

CnipjvAuAip + Cn!iquv*uv£)q- 

In this expression, c,y*/«v are the third order elastic stiff- 
ness coefficients and €,y are the strains resulting from the 
applied stresses. If the applied stresses are within the 
linear elastic regime, the strains are given by 

tij^Sijktvkl (4) 

where s,y*/ are the linear elastic compliances which are 
the inverse of the stiffnesses. 

Thus, if the linear elastic stiffnesses and compliances, 
the third order elastic stiffnesses, and the applied 
stresses are known, an “effective” stiffness tensor can be 
calculated. This can then be used to compute the 
changes in energy flux deviation as a function of applied 
stress which are a result of nonlinear elastic effects. 

The effect of stress on the energy flux deviation was 
modeled for unidirectional T300/5208 gr/ep which was 
assumed to be transversely isotropic. The fiber axis was 
designated to be the X3 axis while the laminate stacking 
direction which is perpendicular to the fibers was 
chosen to be the xi axis. The values of the previously 
measured non-zero independent linear elastic stiffness 
and compliance coefficients for this material are listed 
in Table 1. Likewise, the values of the previously mea- 
sured non-independent non-zero third order stiffness 
coefficient are given in Table 2. In both tables, the con- 
tracted subscript matrix notation is used. 


TABLE 1 


Linear Elastic Stiffness and Compliance Coefficients 


(GPa) 


(GPa)" 1 

Cu 

14.26 

Si 1 

0.092 

Cl2 

6.78 

S 12 

-0.042 

Cl3 

6.5 

S 13 

-0.003 

C33 

108.4 

S 33 

0.0096 

C 44 

5.27 

S 44 

0.190 


TABLE 2 

ThiTd Order (Non-Linear) Elastic Stiffness Coefficients 


Cijk 

(GPa) 

C ijk 

(GPa) 

Cm 

-196 

Cl55 

-49.1 

C 112 

-89 

C 344 

-47 

Ctl3 

-4 

Cl33 

-236 

Cl23 

65 

C 333 

-829 

Cj44 

-33.4 




Calculations were performed for elastic waves propa- 
gating in the X1X3 plane. As in any anisotropic bulk 
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material, three elastic waves will propagate along any 
direction in this plane. Of the three waves propagating 
in this plane, one of them is always a pure mode trans- 
verse (PT) wave with its particle displacement polar- 
ized perpendicular to the xi and X3 axes, i.e., along the 5 
X2 axis. The other two modes are quasi-mode waves 
with components of particle displacements both along 
their direction of propagation and perpendicular to it. 
One is a quasi-transverse (QT) mode wave while the 
other is a quasi-longitudinal (QL) mode wave. All three 1° 
modes suffer energy flux deviation except for propaga- 
tion along the fiber axis (X3) and the laminate stacking 
axis (xi). These are symmetry directions along which all 
three modes are pure mode waves and none suffers 
energy flux deviation. 15 

The energy flux deviation was first computed for the 
condition of no applied stress. Then, calculations were 
performed for two different states of uniaxial stress. The 
first was stress along the fiber direction (X3) of a magni- 
tude of 1 GPa, as shown in FIG. 5. The second was 20 
along the laminate stacking (xj) direction with a magni- 
tude of 0.1 GPa, as shown in FIG. 6. These values are 
near the reported ultimate strengths of this material 
along the respective directions. This allows an estimate ^ 
of the maximum effects of stress on energy flux devia- 
tion. For both conditions, the change in the energy flux 
deviation angle from the condition of zero applied stress 
was computed over the range of propagation directions 
of 0° to 60° from the fiber axis at 2° intervals. 30 

The flux deviation angles of the three modes as a 
function of fiber orientation at zero stress are plotted in 
FIG. 4. The fiber orientation angle is the angle between 
the fibers and the normal N to the wave front or the 
intended direction of propagation. A positive flux devi- 35 
ation angle implies the energy deviates away from the 
fiber direction toward the xj axis while a negative devi- 
ation means that the energy deviates toward the fibers. 
Over this range of fiber orientation angles, the energy of 
the QL and PT mode waves deviates toward the fibers 40 
while that of the QT deviates away from the fibers. 

In FIG. 5, the change in the predicted energy flux 
deviation due to the application of stress along the fibers 
is plotted as a function of fiber orientation angle. The 
energy of the QT mode wave suffers the largest shift in 45 
flux deviation, reaching a maximum of 3° at a propaga- 
tion direction of approximately 20° with respect to the 
fiber direction. The energy of the PT mode wave 
changes by a smaller amount in the opposite direction 
while the QL mode wave suffers a negligible shift. 50 
The relative magnitudes of the flux deviation shifts of 
the different modes can be explained qualitatively by 
considering the ratios of the magnitude of the non-linear 
elastic coefficients to the linear coefficients. The pri- 
mary elastic coefficients affecting the propagation of 55 
the PT and QL modes are those dominated by the ma- 
trix properties. These coefficients are cn, C12, C44, cm, 
and ci 12. The magnitude of the non-linear coefficients 
are over an order of magnitude larger than the linear 
coefficients in this case. However, the ratios of the 60 
non-linear to linear coefficients which dominate the QL 
wave (C33, C133, and C333) are much smaller even though 
the magnitudes of the individual coefficients are larger. 
Therefore, the effect of non-linear elasticity on the en- 
ergy flux deviation should be much smaller for the QL 65 
mode wave. The previous measurements of the effect of 
matrix degradation on energy flux deviation also 
showed a larger change in the flux deviation of the QT 
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mode wave with almost no change in the QL mode 
wave. 

The shift in energy flux deviation due to applied 
stress along the xj axis is shown in FIG. 6. Again the QT 
mode wave suffers the largest change in flux deviation 
angle while the QL mode is almost unchanged. It is 
interesting to note that the direction of the change in 
energy flux is in the opposite direction from the case of 
applied stress along the fiber direction. 

These calculations demonstrate the effect of nonlin- 
ear elasticity on the energy flux deviation of ultrasonic 
waves in gr/ep composite materials. The modes indi- 
cate the angles of fiber orientation and wave modes that 
suffer the maximum shift in flux deviation for the cases 
of applied stress considered. Although the models pres- 
ented were for bulk waves propagating through a thick 
composite material, the same effect is expected for plane 
plate waves propagating in thin plates. The longer prop- 
agation paths possible along plates would make the 
effect more measurable and thus could improve the 
stress resolution. This effect is used to develop a new 
nondestructive method for monitoring stress in compos- 
ite materials or as a new method for measuring their 
nonlinearity. 

In practice, data for flux shifts are generated for a 
specific material having a specific fiber orientation at 
loiown loads. A family of curves for different loads 
could be generated with each curve representing a dif- 
ferent load condition as shown by way of example in 
FIGS. 5 and 6 where only one load condition is de- 
picted. This data and a subsequently determined flux 
shift allows one to approximate the value of an un- 
known load which caused this flux shift. 

Specifically, an ultrasonic wave is applied to a partic- 
ular anisotropic sample material in an initial direction 
and an angle of flux deviation of the ultrasonic wave 
front measured from this initial direction. As discussed 
previously, this flux deviation angle is induced by the 
unknown applied load. Next, the flux shift is determined 
between this flux deviation angle and a previously de- 
termined angle of flux deviation of an ultrasonic wave 
applied to a similar, or the particular, anisotropic refer- 
ence material under an initial known load condition 
such as a no load condition. This determined flux shift is 
compared to a plurality of flux shifts of a similarly 
tested, similar anisotropic reference material, or the 
particular sample material, under a plurality of respec- 
tive, known load conditions, whereby the load applied 
to a particular anisotropic sample material is deter- 
mined. The “s imil ar** anisotropic material(s) are similar 
in that they possess similar fiber orientation, i.e., they 
exhibit flux shifts similar to those of the particular aniso- 
tropic sample material. 

The unknown quantity could likewise be fiber orien- 
tation of a known material, wherein a determined flux 
shift caused by a known load allows for approximation 
of the fiber orientation, as shown in FIGS. 5 and 6. An 
ultrasonic wave is first applied in an initial direction to 
a particular anisotropic sample material under a known 
load and having unknown fiber orientation. Next, an 
angle of flux deviation of the ultrasonic wave front from 
the initial direction is measured. A flux shift is then 
determined between this measured flux deviation angle 
and a previously determined angle of flux deviation of 
an ultrasonic wave applied to the particular anisotropic 
sample material under a known load such as a no load 
condition. This flux shift is then compared to a plurality 
of flux shifts of a plurality of anisotropic reference mate- 
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rials having a plurality of respectively known fiber 
orientations and similarly tested under the same load 
conditions as the particular material, whereby the fiber 
orientation of the particular anisotropic material sample 
is determined. These reference materials are “similar” in 
that they are the same composite as the particular sam- 
ple but have various fiber orientations. 

Many modifications, improvements and substitutions 
will be apparent to the skilled artisan without departing 
from the spirit and scope of the present invention as 
described in the application and defined in the following 
claims. 

What is claimed is: 

1. A method of determining an applied load of a par- 
ticular anisotropic sample material comprising the steps 
of: 

applying an ultrasonic wave to the particular aniso- 
tropic sample material in an initial direction; 

measuring an angle of flux deviation of the ultrasonic 
wave front from the initial direction, the flux devia- 
tion angle induced by the applied load; 

determining a flux shift between this load-induced 
flux deviation angle and a previously determined 
angle of flux deviation of an ultrasonic wave ap- 
plied to a similar anisotropic reference material 
under an initial known load condition; and 

comparing this said flux shift to a plurality of previ- 
ously determined flux shifts of a similarly tested, 
similar anisotropic reference material tested under 
a plurality of respective known load conditions, 
whereby the load applied to the particular aniso- 
tropic sample material is determined by correlating 
said flux shifts to previously determined flux shifts 
due to known loads in said reference material. 

2. The method according to claim 1, wherein said 
measuring step comprises measuring the angle of flux 
deviation in a quasi-transverse wave direction, wherein 
said quasi-transverse wave direction is the direction in 
which particle motion of the wave is predominantly 
perpendicular to the direction of wave propagation. 

3. The method according to claim 1, wherein said 
measuring step comprises measuring the angle of flux 
deviation in a quasi-longitudinal direction, wherein said 
quasi-longitudinal wave direction is the direction in 
which particle motion of the wave is predominantly 
parallel to the direction of wave propagation. 
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4. The method according to claim 1, wherein the 
particular anisotropic material is non-crystalline. 

5. The method according to claim 4, wherein the 
non-crystalline, anisotropic material is a graphite fiber- 

5 reinforced composite. 

6. The method according to claim 5, wherein said 
measuring step comprises measuring the angle of flux 
deviation in a quasi-transverse wave direction, wherein 
said quasi-transverse wave direction is the direction in 

10 which particle motion of the wave is predominantly 
perpendicular to the direction of wave propagation. 

7. The method according to claim 5, wherein said 
measuring step comprises measuring the angle of flux 
deviation in a quasi-longitudinal direction, wherein said 

15 quasi-longitudinal wave direction is the direction in 
which particle motion of the wave is predominantly 
parallel to the direction of wave propagation. 

8. The method according to claim 5, wherein the 
graphite fiber-reinforced composite material is selected 

20 from a group consisting of graphite/epoxy, graphite/- 
magnesium and graphite/aluminum. 

9. The method according to claim 8, wherein said 
measuring step comprises measuring the angle of flux 
deviation in a quasi-transverse wave direction, wherein 

25 said quasi-transverse wave direction is the direction in 
which particle motion of the wave is predominantly 
perpendicular to the direction of wave propagation. 

10. The method according to claim 8, wherein said 
measuring step comprises measuring the angle of flux 

30 deviation in a quasi-longitudinal direction, wherein said 
quasi-longitudinal wave direction is the direction in 
which particle motion of the wave is predominantly 
parallel to the direction of wave propagation. 

11. The method according to claim 1, wherein the 

35 initial known load condition is a no-load condition. 

12. The method according to claim 1, wherein said 
applying step comprises disposing a means for transmit- 
ting an ultrasonic wave on a surface of the anisotropic 
material and said measuring step comprises arranging 

40 means for receiving the transmitted ultrasonic wave on 
another surface of the anisotropic material opposite said 
means for transmitting an ultrasonic wave. 

13. The method according to claim 12, wherein the 
transmitting means is an ultrasonic transducer. 

45 14. The method according to claim 12, wherein the 

receiving means is an array of receiving transducers. 
***** 
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